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ABSTRACT 

High precision cosmology with weak gravitational lensing requires a precise measure of the 
Point Spread Function across the imaging data where the accuracy to which high spatial fre- 
quency variation can be modeled is limited by the stellar number density across the field. We 
analyse dense stellar fields imaged at the Canada-France-Hawaii Telescope to quantify the 
degree of high spatial frequency variation in ground-based imaging Point Spread Functions 
and compare our results to models of atmospheric turbulence. The data shows an anisotropic 
turbulence pattern with an orientation independent of the wind direction and wind speed. We 
find the amplitude of the high spatial frequencies to decrease with increasing exposure time 
as and find a negligibly small atmospheric contribution to the Point Spread Function 

ellipticity variation for exposure times t > 180 seconds. For future surveys analysing shorter 
exposure data, this anisotropic turbulence will need to be taken into account as the ampli- 
tude of the correlated atmospheric distortions becomes comparable to a cosmological lensing 
signal on scales less than ~ 10 arcminutes. This effect could be mitigated, however, by cor- 
relating galaxy shear measured on exposures imaged with a time separation greater than 50 
seconds, for which we find the spatial turbulence patterns to be uncorrected. 

Key words: cosmology: observations - gravitational lensing 



1 INTRODUCTION 

Over the next decade we will see all sky surveys being undertaken 
to address key science questions ranging from cataloging all near 
earth objects to mapping the Dark Univ erse. These surveys include 
the L arge Synoptic S urvey Telescope (|LSST Science Collaboration! 
|200i, LSST), Euclid dUaureiis et al.ll201ll) , and the Panoram i c Sur- 
vey Telescope and Rapid Response System ('Kaise r et"aD l20ld 
PanSTARRS). By observing the growth of dark matter structures 
over cosmic time, weak lensing may provide the most interest- 
ing constraints on the dynamical properties of dark energy. Along- 
side modifications to General Relativity, this dark energy is hy- 
pothesized as an explanation for the accelerated expansion of the 
Universe (see for example Huterer 2010, and references therein). 
The observational measurement of weak gravitational lensing is 
however far from trivial. In fact it is the technique that sets the 
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most stringent requirements on the optical design of future instru- 
mentation as we require accurate knowledge of the Point Spread 
Function (PSF). The PSF also needs to be minimal in terms of 
its size and ellipticity. Dark matter typically induces a one per- 
cent change in the ellipticity of lensed distant galaxies. It is 
this distortion that we measure to constrain cosmology. In con- 
trast the telescope, camera and atmosphere can introduce image 
distortions up to the ten percent level. The main challenge in 
weak lensing is therefore to distinguish the cosmological distortion 
from instrumental and atmospheric distortions in order to recover 
the underlying dark matter signal. A significant effort is under- 
way to improve t he accuracy of weak lensing measurements both 
from a software jKitching et al.l 201Q|) and hardware perspective 
(IlSST Science Collaborationl2009h . Furthermore ambitious space- 
based surveys are proposed in order to remove atmospheric seeing 
and obtain diffraction limited observations dUaureiis et al.ll201ll) . 

For ground-based observatories, the varying refractive index 
in the turbulent layers of the atmosphere produces two effects. It 
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blurs the images of point sources and also introduces an instanta- 
neous ellipticity at the 10% level dde Vries et al. 2007). This can 
be viewed a superposition of speckles. As the exposure time in- 
creases, these speckles de-correlate producing a circular distortion 
pattern. Telescope optics and the camera distortion then induce an 
additional smoothing term and an anisotropic distortion across the 
camera field of view. The combined observed PSF size is roughly 
given by adding the size of the atmospheric PSF in quadrature 
with the size of the instrument PSF. Optimal instrumentation de- 
sign therefore aims to minimize the instrument contribution so that 
the resulting PSF is atmosphere dominated. Image distortions are 
encompassed by the PSF as imaged by stellar point sources. Using 
colour and size information, stars can be separated from galaxies 
and used to model the PSF. For weak lensing measurements to date, 
the PSF model is typically an interpolation between stars, using a 
low order polynomial model to obtain a PSF at each galaxy loca- 
tion. This ass umes however that the PSF do es not vary rapidly on 
small scales (tPaulin-Henriksson et al .112008b . 

In this paper we measure and analyse the image distortion in- 
duced in short exposure observations by atmospheric aberrations. 
This is of particular relevance to the PanSTARRS and upcoming 
LSST ground-based surveys. Both these surveys will image more 
than 20,000 sq. degrees with LSST achieving a depth of r ^ 27.5 
over 10 years. As transient objects are a key science goal for both 
these ambitious surveys repeat short exposure observations are re- 
quired. For LSST each region of sky will be imaged close to 1000 
times with 15 second exposures. The current Pan-STARRS survey 
images each region of sky 5 times with 30 second exposures. As 
both surveys progress a deep data set will be collected for opti- 
mal weak lensing analysis, but with such short exposure times the 
question is whether the rapidly varying atmospheric ellipticity will 
dominate the PSF degrading the ability to model the PSF in indi- 
vidu al ex posures using standard methods. 

IWittman (2005) first looked at this issue by analysing a set of 
10 second and 30 second exposure images of a dense stellar field 
from the Subaru Telescope taken in one night. The conclusion from 
this paper was that the spurious small scale power induced by the 
atmosphere would not be a significant factor in the co-added data 
of future massive surveys based on the finding that the atmospheric 
distortion in exposures separated by 120 s econds was found to be 
uncorrected. This analysis was extended by lAsztalos et al.l (l2007h 
with an analysis of 15 second exposure images from the Gemini 
South Telescope where evidence was presented for a correlation be- 
tween the atmospheric distortion and the ground-based wind speed 
and direction with a consequence that the co-addition of exposures 
taken with the same obs erving site prev ailing wind may not de- 
correlate as concluded bv lWittmanI (l2005h . 

An investigation of the image ellipticity of atmospheric 
aberrations in simu lated short exposure data is presented in 
Ide Vries et al.l (l2007h. Using ray-tracing and Fourier methods they 
simulate a turbulent atmosphere with a single layer Kolmogorov 
power spectrum, finding significant 10% ellipticities in the instanta- 
neous atmospheric PSF. They also investigate the time dependence 
of the atmospherically induced ellipticities. Assuming a constant 
wind speed they show the amplitude of the atmospheric ellipticity 
decreases with time t as t~^/^. This simulation work has been ex- 
tended by Jee & Tyson (2010) with simulated LSST short exposure 
images with a PSF that includes a six layer Kolmogorov power 
spectrum turbulence compon ent and a ray-t r acing optical distor- 
tion model. They recover the Ide Vries et al.l (l2007h result that as 
exposure time increases, the atmospheric PSF becomes rounder. 
Using a Principal Component Analysis (PCA) they model the simu- 



lated short-exposure PSF across the LSST field-of-view and, keep- 
ing only the most significant 20 eigen-modes of that PCA analysis, 
they simulate what is described as a realistic PSF-convolved LSST 
image. They conclude that a typical stellar density is sufficient to 
model the PCA derived PSF and demonstrate their ability to remove 
that PSF distortion from galaxy shape measurement. By limiting 
the PSF to only 20 eigen-modes however, this analysis suppresses 
the random and high spatial frequency variation that we investigate 
in this study. 

In this paper we assess the impact of atmospheric distortions 
on weak lensing surveys by analysing short exposure images of 
high density stellar regions with the wide-field MegaCam Imager 
on the Canada-France-Hawaii Telescope (CFHT). It is of partic- 
ular importance to revisit this issue as a significant improvement 
in the accuracy of weak lensing shape measurement has recently 
been found when data is analysed using individual exposures com- 
pared to the analysis of stacked exposures (Miller et al. in prep). 
Weak lensing data in the future will therefore not be analysed from 
a co-addition of exposures where all atmospheric turbulent effects 
have de-correlated, as it has been to date. With a wide range of 
exposure time, wind speeds and directions we test the findings of 
de Vr ies et al.l (l2007h that the atmospher ic ellipticity decreases with 
time t as t~^/^ and also the findings of lAsztalos et al.l (l2007h that 
the atmospheric ellipticity is correlated with the wind speed and di- 
rection. Our results complement efforts to accurately simulate these 
effects (Chang et al. in prep) and inform the hardware and software 
design of future ground-based surveys accordingly. This paper is 
organised as follows. In Section Owe describe the data and anal- 
ysis methods, presenting the results and analysis in Section [3] and 
conclusions in Section|4l 



2 DATA AND METHODS 

We analyse the CFHT PSF using a series of observations of dense 
stellar fields imaged and archived by the CFHT and listed in Ta- 
ble[T] We select exposures taken in sub-arcsecond seeing conditions 
in the i-band over a range of exposure times. Our main sample of 
60 exposures imaged in 74 seconds were taken during three dif- 
ferent observing runs sampling a wide range of wind speeds and 
directions. We also analyse nine other samples with exposure times 
ranging from 1 second to 450 seconds with each set observed on 
different nights. 

The data were processed and s tellar object catalogues pro- 
duced as described in lHoekstra et al.l (12006) yielding a typical stel- 
lar density of 7 stars per square arcminute, and in all datasets a 
significantly greater stellar density than available in extragalactic 
imaging survey data taken out of the galactic plane. We parameter- 
ize the PSF in terms of stellar ellipticity as measured by weighted 
quadrupole moments 

_ J d'9W{0)i{0)e.ej 

^<pew{0)i{e) ' ^ ^ 

where I{0) is the surface brightness of the star, is the angular 
distance from the star centre and is a Gaussian weight function 
of scale length rg, where Vg is the measured dispersion of the PSF 
dKaiser et al II 19951) . For a perfect ellipse, the weighted quadrupole 
moments are related to the weighted ellipticity parameters by 
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Exptime N Exp Proposal ID Date 



1.0 


28 


08AH22 


2008-04-14/05-07 


2.0 


7 


03BF99 
04BF28 
03BF09 


2004-01-19/23 
2004-12-05/12 
2003-11-18/12-16 


7.5 


3 


08AQ98 


2008-05-09/06-03 


10.0 


33 


06BH48 


2006-08-19/22 


30.0 


9 


03BF99 


2004-01-19/23 


45.0 


5 


03BQ97 


2003-09-30 


74.0 


60 


08AH22 


2008-04-14/05-08 


180.0 


9 


03BL06 


2003-10-03 


300.0 


11 


08AH53 


2008-02-13/14 


450.0 


12 


07BC02 


2007-11-14/01-10 



Table 1. Table of CFHT archive data used in this analysis listing the ex- 
posure time, total number of exposures used, the CFHT Proposal ID which 
can be used to query the CFHT Science Data Archive, and the range of 
observation dates. 



where R is related to object size and given by 

R = VQii + Q22 , (3) 

If the weight function W{0) = 1 in equation [T] the ellipticity or 
polarisation \s\ = ( 1 — + 0^), where (3 i s the axial ra- 

tio of the ellipse fsee lBartelmann & SchneideJl200lh . We note that 
for the purposes of this study the quadrupole moment description 
of the PSF that we adopt is sufficient. For a detailed weak lensing 
analysis however, many techniques now determine either a set of 
orthonormal 2D basis functions or a 2D pixel-based model to de- 
scribe the PSF ( Kitching et al. 2010). Both these methods would 
be more challenging to model in the presence of high order spatial 
frequency variation. 

In each exposure the variation of the PSF ellipticity sf^^ and 
size R across the field of view is considered to have two compo- 
nents; a smoothly varying second order polynomial over position, 
^modeij^^^ ^) and R{x, y) for each chip and higher spatial varying 
residuals with 



OSi — Si 



SR^ — RpsF 



td2 

-'^model ■ 



(4) 



(5) 



With a typical number density of stellar objects imaged in an high- 
Galactic latitude field (40 per MegaCam chip, compared with over 
700 in this analysis), the chip- wise second order polynomial model 
would be the most complex model that could be accurately fit to 
lensing survey data CRowe.2010.) . 



3 ANALYSIS AND RESULTS 

Figure [T] shows a typical PSF pattern for a 74 second exposure as 
imaged by MegaCam at CFHT. The field of view is a square degree 
with a pixel scale of 0.186 arcseconds and an 9 x 4 CCD pattern 
which is visible. The left and right panels show the variation in the 
amplitude of the average PSF si and S2 respectively across the field 
of view. Each grid point in the greyscale map contains on average 5 
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Figure 1. A typical PSF pattern for a 74 second exposure. The left and right 
panels show the variation in the amplitude of the average PSF £i and £2 
respectively across the field of view. The upper panels show the observed 
ellipticity variation, the middle panels show the second order polynomial 
model fit to the data. The lower panels reveal the residual ellipticities dei 
and 6 £2 components. The residuals show high spatial frequencies and a 
preferred direction. The ellipticity amplitude in each panel is indicated by 
the greyscale shown in the vertical colour bar. 



stars and spans 0.3x0.3 arcmins. In an image of a typical field, out 
of the Galactic plane, we would find ^ 1 useable star in every ^ 5 
grid points (^0.4 stars per sq arcmin). The upper panels show the 
observed ellipticity variation. The middle panels show the second 
order chip- wise polynomial model fit to the data. The lower panels 
reveal the residual ellipticities Sei and Se2 components. Figure |2] 
shows the size variation for the same exposure with the size resid- 
uals (right hand panel) showing the same high spatial frequencies 
and a preferred direction as the ellipticity residuals. For the rest 
of the paper we focus mainly on the ellipticity variation, but come 
back to the size variation in Section|4]where we show the high spa- 
tial frequency variation of PSF size is as detrimental to the weak 
lensing shape measurement as the variation in ellipticity. 

In order to investigate the high spatial frequency variation of 
the PSF we measure the two-point correlation function of the resid- 
ual PSF ellipticities in Figure [3] showing the average systematic 
residual PSF correlation functions ^+ (upper two panels) and ^_ 
(lower panel) where 



9')et{e' + e)) ± {er{e')sr{e' + e)) (6) 
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5000 10^1.5x10^ 




1.5x10^ 



Figure 2. A typical PSF size pattern for the same 74 second exposure shown 
in Figure [T] The left hand panel shows the variation in PSF size R across 
the field with the grey scale in arcseconds. The right hand panel shows the 
residual variation in size after a second order polynomial model has been 
removed. The ripple pattern in size follows the same structure as the ellip- 
ticity residuals. The right hand grey scale is in units of 10~^ arcsec^ . 



and st,r are the tangential and rotated ellipticity parameters rotated 
into the reference frame joining each pair of correlated objects. The 
average is taken over all exposures in our sample split by exposure 
length. Figure [3] shows the amplitude of the residuals for four sets 
of decreasing exposure time revealing a characteristic shape to the 
correlation function that we find for all exposure times. <J+ is pos- 
itively correlated < 2\ anti-correlated 2' < < 7\ then pos- 
itively correlated until it becomes consistent with zero > 15'. 
This c haracteristic shape is also seen in the results of IWittmanI 
(l2005h . ^_ is measured with significance only for the shortest expo- 
sures (shown lower panel) and is positively correlated for all scales 

< l5^ 

The residual PSF correlation functions measured between 
consecutive Is, 10s and 74s exposures were found to be consistent 
with zero. We can therefore conclude that the atmospheric turbu- 
lence distortion de-correlates in < 50 seconds where the timescale 
is set by the CFHT MegaCam read-out, o verheads, sle w and ac- 
quisition time. This is in agreement with WittmanI (|2005|) who set 
a de-correlation time < 120 seconds, limited by the Subaru read 
time. 



3.1 Comparison to a simple, isotropic turbulence model 

The von Karman model for isotropic atmospheric turbulence pre- 
dicts a projected, two-dimensional power spectrum 



P{1) oc 



+ ■ 



-11/6 



(7) 



whe re the angle Oq defines the outer scale of the spectrum (e.g. 
ISasi ela 1994). In this model, a physically-motivated modification 
of the scale-free Kolmogorov spectrum, turbulent structures in im- 
ages of the sky become uncorrected at separations greater than Oo . 
As shown in Appendix lAl the isotropic ^+ correlation function for 
such a spectrum may then be written as 



e+(^)oc^'/'i^_5/6 (27r^/^o), 



(8) 



where Ku(x) is the modified Bessel function of the second kind 
lArfken & Weber 2005). Numerical approximatio n of thes e func- 
tions is simple via the useful series expansion of iKostrounI Ji98Q|) . 
The functional form of ^- is given in equation (IA4I) . This expres- 
sion is significantly more complicated than the ^+ case and was 
found to be prone to numerical instabilities. 

We now investigate whether the model given in equation JS]) 




+ 
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[arcmin] 

Figure 3. The ^+ {6) and ^_ {6) correlation function estimates for the PSF 
ellipticity data, with best-fitting von Karman models overlaid on the ^+ re- 
sults and the corresponding ^_ prediction plotted over the measured ^_ 
results (see Section [TTI and Appendix |A). Upper panel: the t = Is (solid 
line: best-fitting model) and t = 7.5s (dashed line: best-fitting model) 
^+{0) results. Middle panel: the t = 74s (solid line: best-fitting model) 
and t = 450s (dashed line: best-fitting model) ^+{0) results. Lower panel: 
the t = Is, and t = 7.5s, ^_ (9) results (the longer exposure PSF patterns 
showed ^_ consistent with zero on all scales). 



is able to reproduce the correlation function results seen in Fig- 
ure [3] We restrict ourselves to fitting only the ^+ model, which 
can be simply calculated and carries almost all the detectable sig- 
nal. We go on to compare predictions for ^_ based on these fitting 
results to measurements of ^_ in the data. In Figure [3] we have 
overlaid the best-fitting von Karman models to the ^+ correlation 
function estimates for each of the example t — Is, 7.5s, 74s & 
450s exposure da ta sets. The Levenberg-Marquardt algorithm (e.g. 
Press et al was used to determine the maximum-likelihood 

parameter fits to the overall amplitude of the expression in equation 
([U, and for the outer scale length Oq. In the lower panel of Figure[3] 
we use the best-fitting values of Oq from the ^+ models to gener- 
ate a von Karman prediction for the ^_ signal in the t=ls and 7.5s 
exposure data. 

It can be seen that the von Karman model gives a reasonable 
fit to the inner slope of ^+ for the three t = Is, 7.5s & 74s ex- 
posure data sets in Figure [3] For these three sets we also found 
broadly consistent best-fitting values of the outer scale in the range 
Oq — 2.62-3.22^ These results suggest a turbulent origin for the 
small-scale correlated ellipticity measurements in these image, but 
in all three cases the von Karman model fails to capture both the 
trough on scales ^ ^ 4^ and the second peak in correlation near 
10^ The discrepancy is also marked in the comparison between the 
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^_ model and the data. Here a very distinct excess 'bump' is seen, 
again peaking around ^ ~ 4^ in the Is exposure time data. How- 
ever, as the ampUtude of the ^- correlation is consistently smaller 
than that of these results are noisier. We discuss potential ori- 
gins for these features in Section]?] For the t = 450s data we found 
a model amplitude consistent with zero, indicative of the weakness 
of the signal for this long exposure data. 



3.2 Dependence of the amplitude of the atmospheric 
distortion on exposure time 

A key result from the simulations of Ide Vries et al.l (l2007h showed 
that the ellipticity introduced by atmospheric aberrations decreases 
with time as t~^/^ assuming a constant wind speed. The two-point 
correlation function of the atmospheric variations ^+ is therefore 
expected to decrease in amplitude as t~^. Figured shows the aver- 
age two-point correlation function |^+| as a function of the expo- 
sure time. Results are shown for two representative angular scales 
with ^+ measured at ^ = 0.7^ (upper) where the von Karman at- 
mospheric turbulence model provides a good fit to the data, and at 
= 4.9' which corresponds to the first negative dip in the correla- 
tion function. We find that the de Vries et al. ( 2007) model (shown 
dot-dashed) is a good fit to the data but note that for this analysis, 
we were unable to select a significant sized sampl e of exposures 
with a constant wind speed as is assumed in the Ide Vries et al.l 
theory. We conclude that the scatter we see in the results 
must be in part driven by the range of wind speeds within our sam- 
ple. 

Figure |4l also allows us to compare the amplitude of the at- 
mospheric turbulence to the expected cosmological signal in the 
CFHT Lensing Survey (CFHTLenSQ). For both angular scales 
shown we over-plot two horizontal bars that indicate the ampli- 
tude of the two-point shear correlation function ^+{0) for two low 
redshift tomographic bins with a mean reds hift of z — 0.37 an d 
z — 0.54 assuming a WMAP7 cosmology jKomatsu et al.ll20ll]) . 
We should note here that one cannot directly relate the PSF ellip- 
ticity correlation with that expected from shear, as the effect of 
the PSF on the shear measurement depends on the relative size 
of galaxy such that it would negligible for large galaxies or po- 
tenti ally amplified for galaxie s comparable to or smaller than the 
PSF dPaulin-Henriksson et al.ll20b 8). This comparison does how- 
ever demonstrate that even for the lowest redshift tomographic bin, 
the residual systematic error from atmospheric turbulence is well 
below the cosmological signal for the 600 second CFHTLenS ex- 
posures. 



3.3 Dependence on wind speed and direction 

In order to compare the turbulence patterns observed to the wind 
direction we quantify the direction of the atmospheric turbulence 
using correlation as a function of vector separation and direction 
on the sky, rather than the commonly-used, azimuthally- averaged 
functions ^+{0) = (J+(|0|). We calculate the ^+ residual correla- 
tion as a function of vector angular separation 9 using the Discrete 
Fourier Transform to find the power spectrum of the 5ei and 5e2 
images (examples of which are shown in Figure[T]), then employing 
the Wiener- Khinchin theorem to generate images of ^+{0). 



^ CFHTLenS analyses data from the CFHT Legacy Survey, using both the 
Wide, Deep and Pre-imaging surveys 




1000 



exposure time (sees) 



Figure 4. The evolution with exposure time of the atmospheric turbulence 
as measured by the two-point correlation function ^+ at two different angu- 
lar scales; 0.7 arcmin(black upper) and 4.9 arcmin (red lower). The best-fit 
de Vries et al.(2007) prediction is shown dot-dashed for each angular scale. 
This can be compared to the two pairs of thick horizontal bars over plot for 
each angular scale. These bars indicate the expected CFHTLenS cosmolog- 
ical signal for two tomographic bins with a mean redshift z = 0.37 and 
z = 0.54. determining the exposure time below which the cosmological 
signal becomes lower than the residual atmospheric PSF signal. 





-20 20 

6^ (arcmin) 



-20 20 

(arcmin) 




-20 20 

6^ (arcmin) 



-20 20 

6^ (arcmin) 



Figure 5. The ^+ residual correlation as a function of 2D vector separation 
6 for four example 74 second exposures. The dominant direction of the 
PSF residuals can be compared to the wind direction shown as an arrow 
inset. The upper left panel shows an example exposure with a wind speed of 
14ms~^. The upper right panel shows the exposure also shown in Figure[T] 
and has a wind speed of 6ms The lower left panel has a wind speed of 
2ms ~ ^ and the lower right panel has a wind speed of llms~ ^ . 



Four example results are shown in Figure \5\ revealing a non- 
isotropic function. The example exposure shown in Figure [T] is the 
same exposure presented in the upper right hand panel of Figure [5] 
A comparison of these two Figures shows that the elliptical ^+ ( 1 1 ) 
has a central dipole with a major axis along the direction of the rip- 
ples seen in the 6ei maps of Figure [T] For other exposures however 
we see what appears to be a superposition of ripple patterns (see for 
example the upper left and lower right panels of Figure [5]), poten- 
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Figure 6. A comparison of wind direction and PSF residual direction as a 
function of the ground wind speed measured at the start of each observation 
for three samples grouped by exposure time with t ^ 10s (stars), 10s> t ^ 
45s (circles) and t = 74s (crosses). 

tially arising from different turbulent layers in the atmosphere. In 
these cases the orientation of the central dipole indicates the aver- 
age ripple direction. Based on these comparisons over the full data 
set we use the orientation of the central dipole to determine an ef- 
fective ripple orientation and we measure this orientation using the 
quadrupole moments in equation [3 with rg = 2.3 arcmin. Apply- 
ing this method to all exposures with t ^ 74s, where the amplitude 
of the atmospheric turbulence is sufficiently high, we can then com- 
pare this ripple direction to the wind direction relative to the image 
(shown by the arrow in the upper inset in Figure [5]), calculated as 
described in Appendix iBl 

Figure [6l shows a compilation of results for three sets of data; 
t ^ 10s (stars), 10s> t ^ 45s (circles) and t = 74s (crosses) 
comparing wind direction and PSF or ripple orientation. We find 
that in general there is little evidence of strong correlation between 
wind direction and PSF orientation for the r ange of exposure times 
tested. This result is in disagreement with lAsztalos et al.l (l2007l) 
who found a relationship between wind direction and PSF ellip- 
ticities for wind speeds ranging from 2-6 m/s over 4 consecutive 
nights. We note that we could have drawn a similar conclusion had 
we analysed a smaller set of exposures, as the points that cluster 
in Figure [6l are typically taken on the same night. This clustering 
with observation date, if sets of observations were to be analysed 
in isolation, might well suggest a preferred wind-PSF orientation 
which we do not find when we analyse a large set of data spanning 
many years of CFHT imaging. 



4 CONCLUSIONS: IMPACT FOR WEAK LENSING 
OBSERVATIONS 

In this study we have analysed short exposure observations of dense 
stellar fields to quantify the high spatial frequency variation of the 
PSF which we attribute to atmospheric effects. We have shown that 
in short exposures with t ^ 30s, the atmosphere contributes a sig- 
nificant correlated anisotropy to the PSF on angular scales ^ < 10^ 
that would dominate the cosmological signal measured in the lower 
redshift bins of a cosmological tomographic analysis. On these an- 
gular scales the high spatial frequency of the atmospheric aberra- 
tion is too rapid to model with a typical stellar density and standard 
methods. 




10 



theta (arcmin) 

Figure 7. The measured residual ellipticity correlation in our t = 10s data 
set compared to the expected WMAP7 cosmological two-point shear cor- 
relation signal for three tomographic bins with a mean redshift z = 0.37, 
z = 0.54 and z = 0.79 (solid lines, the lowest amplitude corresponds to 
the lowest redshift bin). The data oscillates between positive and negative 
correlation which we indicate on this log-log plot using filled points where 
the data is negative which is primarily between the two dashed vertical lines. 



This does not mean, however, that multiple short exposure im- 
ages cannot not be combined in a way that reduces the cumula- 
tive impact of this (random) atmospheric anisotropy on estimates 
of gravitational shear in the final analysis. It was found in this 
study that the turbulent patterns in consecutive exposures are sta- 
tistically uncorrected when separated by timescales ^ 50s. Our 
observations motivate the combination of information from mul- 
tiple short exposure images in a way that takes advantage of this 
fact. Examples of such approaches might be cross-correlating shear 
estimates from different exposures to estimate the shear autocor- 
relation, or combining images into a stacked, co- added image. If 
stacking, great care must taken to retain control and knowledge of 
chang es in the PSF due to interpol ation and stacking of dithered 
data. Rowe, Hirata & Rhodes 2011 present an example of a linear 
image stacking approach that seeks to preserve this information. 

As well as finding that consecuti ve exposures are unco rrelated, 
we have confirmed the predictions of Ide Vries et al.! (l2007h in find- 
ing that the amplitude of atmospheric distortion patterns decreases 
with time as t~^/^. Our results also show that this effect is not a 
significant source of error for the CFHTLenS survey where for the 
first time a lensing survey is analysing the 600 second exposures in- 
dividually (M iller et al in pr ep), as compared to previous analyses 
using a stack (IFu et al.ll2008h . 

Figure[7]compares the measured residual ellipticity correlation 
in our t = 10s data set with the expected cosmological signal for 
three tomographic bins with a mean redshift z = 0.37, z = 0.54, 
z = 0.79. This is the closest data set analysed to the proposed 
t = 15s exposures for LSST. This figure shows that our findings 
agree to some extent with a Jee & TvsonI (I201G) , who focus on an- 
gular scales ^ > 10^ and find that the PSF can be modeled to high 
accuracy on these scales. For smaller angular scales, however atmo- 
spheric turbulence effects will be a significant source of systematic 
error in a weak lensing analysis if not corrected for. Furthermore, 
as a result of the anisotropic nature of the residual correlation func- 
tion (shown in Figure [5]) the effect of atmospheric turbulence may 
leak to larger scales. 

For a futu re generation all- sky lensing survey 

lAmara & RefregieJ ilOOli) derive a requirement on the vari- 
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ance of systematic errors to be below a^ys < 10 ~^ such that 
experiments are Umited by statistical noise rather than systematic 
errors (see also van Wa erbeke et al. (2006) for a simil ar conclu- 
sion). To reach this goal IPaulin-Henriks son et alJ (l2008h show that 
this requires cr [spsf] ^ 10~^ for each ellipticity component and 
c^[^psf] ^ 10"^i?i>sF for ^PSF ^ 0.05 and a typical galaxy/PSF 
size ratio of Rgai/RpsF ^ 1.5. For the t = 10s data we measure 

^[^psf] = = (9.9 ± 1.7) x 10"' (9) 

^i^ = ^|f^ = (9.6±2.0)xl0- (10) 

where the error is calculated from the scatter between the 33 ex- 
posures in the t = 10s data set showing the variance within the 
set, not the error on the mean. These measurements are an order 
of magnitude larger i n both size and ellip t icity c ompared to the re- 
quirements quoted in lAmara & Refregiei^ (lOOH). 

Using equation (15) inlPaulin-Henrikss on et al ] (l2008h we can 
calculate the variance of the systematic errors we would expect 
from the residual atmospheric size and ellipticity variation that we 
measure from the t = 10s data set: we find cTgys ^ 10~^. This 
value can be interpreted as the systematic variance produced in 
multiple 10s exposures if: (i) the turbulent PSF cannot be mod- 
elled due to insufficient stellar density; and (ii) no attempt is made 
to combine shear estimates from multiple exposures for which the 
turbulence patterns are mutually statistically independent. We note 
that this calculation depends on the accuracy of both PSF elliptic- 
ity and PSF size characterisation, and both these are affected by 
atmospheric turbulence. 

We note that while the expectation of the turbulent contribu- 
tion to £psF across multiple exposures is zero, assuming any static 
component is removed, this is not the case for i?psF- This has im- 
plications for plans to cross-correlate using only shear estimates 
from different exposures to suppress PSF model errors. The utility 
of such a scheme is that the ellipticity cross terms {SsiSej) (for 
i / j) have zero expectation value, and the variance of this prod- 
uct is the square of the individual variances. This product variance 
decreases oc l/A/gxp, where A/exp is the total number of exposures 
being cross-correlated, and therefore cr[spsF] decreases rapidly as 
oc 1/A/exp in the cross-correlation. 

However, the same approach does not reduce uncertainty in 
RpsF same rate because this quantity has an unknown ex- 

pectation value which must be characterized. Instead, cr[i?psF] ^ 
1 / i/A/exp. This differing behaviour in the combination of multiple 
exposures must be considered when comparing the results in equa- 
tions ^ & ([TOl ) to the required levels for the estimation of cosmic 
shear. The characterization of turbulent variation in Rps^ there- 
fore assumes added importance for upcoming surveys. Combined 
with A/exp, the single-exposure value quoted in equation ( [TOl ) can 
be used to estimate this important effect for surveys using multi- 
ple exposures of duration t ^ 10s, assuming similar atmospheric 
conditions to those at CFHT. 

Both the results in this paper and lWittmanI (l2005h show an os- 
cillatory pattern in the residual ^+ correlation function at ^ > 2.5 
arcminutes, a pattern which is not consistent with the model pre- 
dicted by a simp le, isotropic von Karman spectrum (see Figure [3]) . 
IWittmanI (l2005h describes this ringing as an artifact of the inter- 
polation scheme used to model the PSF on each individual chip 
which for the Suprime-Cam camera used in these observations has 
a chip size of 13.65^ x 6.83^ In this analysis we apply a similar, 
low-order chipwise PSF model, where for MegaCam the chip size 



is 12.70^ X 6.35^ We thus agree with lWittmanl ( l2005h that one pos- 
sible origin for the 'ringing' is over- subtraction of a chipwise PSF 
model for the telescope optics, as the first dip occurs at a close to 
half the shorter dimension of the chip. 

However, it is less clear why chipwise fitting would cause an 
excess at 10^ In addition, looking at the ^+(9) (i.e. not circularly 
averaged) two-point correlation functions (e.g. Figure [5]) we see 
that the oscillatory pattern is anisotropic in nature but often not 
aligned with the x-y chip grid. This alignment would be something 
to be expected if it were an artifact of the chipwise PSF model re- 
moval. We therefore offer an alternative suggestion that these rip- 
ples have an atmospheric origin, related to anisotropy in the atmo- 
spheric power spectrum (c.f. Figure [5]), which cannot be generated 
by simple isotropic models such as the von Karman. This con- 
clusion is supported by the finding that the residual patterns are 
uncorrected between exposures. However, without further data, or 
access to imaging with chips of differing dimensions, it remains 
difficult to say definitively which interpretation is correct. 

Although not performed explicitly in this study, it is worth- 
while mentioning the separation of the correlated turbulence sig- 
nal into 'E' and 'B' mode contributions (see, e.g., Crittenden et al.l 
l2002h . The atmospheric PSF pattern for a turbulent sky may con- 
tain both E- and B -modes; if this PSF is not properly modelled 
any shear catalogue will then contain an imprint of this pattern. An 
E/B-mode decomposition of the cosmic shear signal on the scales 
where atmospheric turbulence patterns are correlated might there- 
fore detect this effect in the form of a non-zero B-mode. However, 
telescope optical abberrations, which also contribute strongly to 
PSF anisotropy across most angular scales of interest fo r lensing, 
tend to generate a more significant E-mode than B-mode jHoekstral 
l2004l) . In addition, cosmological effects can also generate B -modes 
fe.g. JSchneider et al.l2002l : IVale et al .120041) . The absence of strong 
B -modes in cosmic shear measurements therefore remains a prob- 
lematic indicator of the absence of systematic errors due to PSF 
anisotropy, and the potential presence of B -modes in turbulent PSF 
patterns does not alter this fact. 

Finally, when considering our full sample of data, we have 
found little evidence for a strong and consistent relationship be- 
tween the preferred direction of anisotropic correlations in atmo- 
spheric distortion and the ground-based wind velocity. We do not 
necessarily expect to find such a relationship, as the all the tur- 
bulent layers responsible for the atmospheric PSF pattern are not 
necessarily correlated with the ground conditions. However, we do 
find that consecutive observations seem to yield similar angular off- 
sets between wind direction and turbulent PSF pattern orientation: 
such effects might lead to an over-interpretation of the correlation 
between wind direction and PSF pattern orientation in more limited 
data sets. In all, these results do not strongly motivate the use of the 
ground-based wind direction and speed as an input to a PCA-type 
model of the atmospheric part of the PSF on short exposure data. 
However, it may still be useful in modelling out repea table PSF dis- 
tortions due to telescope motion under wind (see, e.g. Jjarvis & JainI 
2004). Once again, this motivates data analysis strategies that take 
advantage of the turbulent pattern being uncorrected between suc- 
cessive exposures. 
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APPENDIX A: CORRELATION FUNCTIONS FOR THE 
VON KARMAN POWER SPECTRUM 

As in the case of cosmological s hear, the ellipticity correlation 
functions ^+(6*) and ^-(0) fsee IBartelmann & Schneidej l200lh 
due to atmospheric effects can be related to power spectrum of the 
von Karman model as 



^+(6*) oc / ldlJo{27TW) 
Jo 

^-{0) oc / mjA{2nlO) 
Jo 



-11/6 



-11/6 



(Al) 
(A2) 



where Ju{x) is the i/-th order Bessel function, and where we are 
assuming that the spatial power spectrum of ellipticity distortions 
due to turbulence follows P(/) as given in equation 0. Using these 
expressions, we calculate the ^+ correlation function as 



(A3) 



where K^{x) is the mo dified Bessel function of the second kind 
(lArfken& WebeiJl2005h . The expression for ^_ is more compli- 
cated: 



^-7/6 

If, 

^-7/6 

IStt 



9§ + 7rV)7_i7/6(27re/eo) 



-1/6 J 
-1/6, 



-ll/i 



^23/6 



385 
432 



Ol- 707T^O^ r(-7/6) 



(A4) 



where Iu(x) is the modified Bessel function of the first kind and 
T{x) is the Gamma function ( Arfken & Weber 2005). The func- 
tions Iu{x) diverge rapidly for x > 1, and the expression in equa- 
tion (IA4b shows signs of numerical instability for > 40o. In this 
study we therefore make a crude approximation to avoid this be- 
haviour. We manually set ^_ (0) = where > 4.20o for the von 
Karman power spectrum model when investigating or plotting ^_ . 



APPENDIX B: THE PROJECTED DIRECTION OF 
GROUND WIND IN IMAGES ON THE CELESTIAL 
SPHERE 

We consider two general points with altitude- azimuth coordinates 
(ai, Al) and (a2, A2) in the horizontal coordinate system, which 
uses the observer's local horizon as the fundamental plane. As for 
any two points on the coordinate sphere, a single great circle may be 
constructed that passes through (ai, Ai) and (a2, A2). We label as 
a the angular size of the arc between these two points on the great 
circle. Using the spherical law of cosines, it is straightforward to 
show that a is given by 

cos (cr) — COS (ai) COS (a2) COS (Al — A2) 

+ sin(Ai)sin(A2). (Bl) 

This result will be used below to derive the projected direction of 
ground wind on telescope images. 

The pixel ^-axis of the MegaCam instrument mounted on 
CFHT rotates to align with the celestial North Pole. In the hor- 
izontal coordinate system defined at the telescope, celestial North 
has the altitude- azimuth coordinates N = (0,0) where is the geo- 
graphic latitude of Mauna Kea. We define the telescope pointing di- 
rection (i.e. the direction of the image field centre) as T = (at, At) 
in this coordinate system. The points N and T on the celestial 
sphere can be connected by the arc of a great circle of constant 
Right Ascension. The y-3xis of CFHT images will align with this 
Meridian. 

The direction of ground wind is measured at CFHT and sup- 
plied with the astronomical images in degrees relative to compass 
North along the ground. In the horizontal coordinate system at 
CFHT this direction can therefore be expressed as W = (0, Aw), 
where Aw is the measured wind direction. We may then define a 
second great circle passing through W and T. The angle at which 
this great circle intersects the Meridian through T and N is there- 
fore the projected direction of ground wind on the telescope image, 
relative to the CCD ^-axis. 

The spherical law of cosines gives the following result for this 
projected wind direction angle 0: 



cos (0) 



cos (7) — cos (a) cos (/3) 



sin (a) sin {f3) 
where from equation dBlb we have 



(B2) 



cos (a) = cos (at) cos (Aw — At), (B3) 
cos (/3) = cos (at) cos ((/)) cos (At) + sin (at) sin (0), (B4) 
cos{'y) — cos (Aw) cos (0). (B5) 

However, the expressions above do not provide the sign of 0, i.e. 
whether W lies to the West or East of the Meridian through T and 
N, the 2/-axis on the telescope image. This information is needed 
to correctly calculate the angular separation between the wind di- 
rection and primary PSF residual direction as described in Section 
[331 

The positive direction on the CCD x-axis of MegaCam can 
be associated with positive hour angle H in the equatorial coordi- 
nate system. If we define the positive direction as that moving 
clockwise from the y-axis, then it is straightforwardly seen that 



sgn(^) = sgn(ifw - Ht) 



(B6) 



where sgn(x) is the signum function, is the hour angle of the 
wind direction vector W and Ht is the hour angle of the telescope 
pointing direction T at the local sidereal time of observation. The 
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telescope hour angle Ht is provided by CFHT as metadata along- 
side all astronomical images, so it remains solely to determine H^. 

Using standard formulae for the conversion of horizontal to 
equatorial spherical coordinates, the declination of the wind vec- 
tor W is given by 



sin (^w) = cos {(/)) cos (Aw) 
The hour angle of W is then given by 



(B7) 



sin (Aw) 



cos (5w 



— arcsm 



sin (Aw 



V^l -cos2 (0) cos2 (Aw) 



(B8) 



When substituted into equation (IB6I) this allows the sign of to be 
calculated. In combination with the determination of the magnitude 
of from equation (IB2I) , the projected direction of ground wind on 
a telescope image pointing in the direction T is fully specified. 
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